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ABSTRACT 



Aims. We investigate the atmospheric structure and fundamental properties of the red supergiant VY CMa. 

Methods. We obtained near-infrared spectro-interferometric observations of VY CMa with spectral resolutions of 35 and 1500 using 
the AMBER instrument at the VLTI. 

Results. The visibility data indicate the presence of molecular layers of water vapor and CO in the extended atmosphere with an 
asymmetric morphology. The uniform disk diameter in the water band around 2.0 //m is increased by ~20% compared to the near- 
continuum bandpass at 2.20-2.25 yum, and in the CO band at 2.3-2.5 /jm it is increased by up to ~50%. The closure phases indicate 
relatively small deviations from point symmetry close to the photospheric layer, and stronger deviations in the extended H2O and 
CO layers. Making use of the high spatial and spectral resolution, a near-continuum bandpass can be isolated from contamination 
by molecular and dusty layers, and the Rosseland-mean photospheric angular diameter is estimated to 11.3±0.3mas based on a 
PHOENIX atmosphere model. Together with recent high-precision estimates of the distance and spectro-photometry, this estimate 
corresponds to a radius of 1420± 120 R and an effective temperature of 3490+90 K. 

Conclusions. VY CMa exhibits asymmetric, possibly clumpy, atmospheric layers of H 2 and CO, which are not co-spatial, within a 
larger elongated dusty envelope. Our revised fundamental parameters put VY CMa close to the Hayashi limit of recent evolutionary 
tracks of initial mass 25 M Q with rotation or 32 M Q without rotation, shortly before evolving blueward in the HR-diagram. 
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1. Introduction 

VY Canis Majoris (VY CMa) is one of the most luminous and 
most massive red supergiants, and thus most likely progenitors 
of a core-collapse supernova (SN) in our Galaxy. Because of its 
brightness and importance, it has traditionally been extensively 
observed at many wavelengths. Currently, VY CMa experiences 
a renewed and increased interest because of observations at new 
facilities in cluding, for instance, the Herschel sp ace observator y 
dRover et a l. 2010) and the submillimeter array (iFu etal.1 12012). 
as well as because of new results regarding the complexity of its 
molecular envelope (e.g., Ziurys et al. 2 0071) . 

While is has been clear that VY CMa's fundamental pa- 
rameters are extreme, their exact values have been debated dur- 
ing the last decade with luminosities between 6 x 10 4 L Q and 
5 x 10 5 L Q , effective temperatures between 2700K and 3650K, 
radii between 600 R a and 3000 R a , and initial masses between 
12 M and 40 Mq dSmith et alJ200lt iMonnier et all2000L l2004t 
iMassev et al.l l2006l) . Very rec ently, precise paral lax measure- 
ments bv lchoi et all d2008l) and lZhang et alj d2012l) significantly 
improved the distance estimate to about 1 .2 kpc, which together 
with r ecent accurate broad-band photometry by Smith et al.l 
(2001) points to a luminosity of about 3 x 10 5 L G . 

VY CMa is embedded in a large optical nebula (Herbig 
1 1972b and is al so an extended object in the infrared 



dMcCarthv & Low! 1 19751) . Several studies have investigated the 



* Based on observations made with the VLT Interferometer (VLTI) 
at Paranal Observatory under programme ID 386.D-0012 



asymmetric circumstellar environment in detail and revealed a 
complex asymmetric morphology (e.g. Wittkowski et al. 1998; 
Monnier et al. 1999, 2000, 2004; Smith et al. 2001; Muller et al. 
2007; Humph reys et all I2OO7P iJones et al.l 120071) . Recently. 
Smith et al. (2009) proposed a scenario of VY CMa's circum- 
stellar environment that consists of fast and dense CO/dust 
cloudlets that are ejected into a larger relatively slow and less 
dense asymmetric envelope due to a strong and variable mass 
loss. 



2. Observations and data reduction 

We observed VY CMa with the VLTI/AMBER instrument in 
visitor mode on 6 and 7 March 2011. Table Q] lists the details 
of our observations. We used the E0-G0-H0 configuration of 
the VLTI, providing projected baseline lengths between 11.2m 
and 48.0m. We obtained 5 observations using the low reso- 
lution mode with a spectral resolution of R ~ 30 at differ- 
ent hour angles, thus different projected baseline lengths and 
angles. We also obtained one observation using each of the 
two /f-band medium resolution modes, together covering wave- 
lengths between 1.92/im and 2.47 //m at a spectral resolution of 
R ~ 1500. All VY CMa observations are embedded between 
observations of an interferometric calibration star. We chose 
HR 3052 (spectral type K5 III, angular diameter 2.62+0. 19 mas) 
and 1 Pup (K5III, 3.85+0.28 mas) fr om the ESO C alibrator 
Selector CalVin based on the catalog bv lLafrasse et al.l (l2010l) . 
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Table 1. Log of our VLTI/ AMBER obsei-vations 



Target 



Date Time 



Mode 



2011-03 



UT 



Baseline PA 
m dee 



VYCMa 
VYCMa 
VYCMa 
VYCMa 
VYCMa 
VYCMa 
VYCMa 



-06 

-06 
-06 

-07 
-07 
-07 
-07 



1:42 
3:05 
4:14 
0:55 
1:35 
3:02 
4:29 



Low-JHK 
MR-K 2.3 
Low-JHK 
Low-JHK 
Low-JHK 
Low-JHK 
MR-K 2.1 



15.8/31.6/47.5 
14.4/28.7/43.0 
11.8/23.6/35.4 
16.0/32.0/48.0 
15.9/31.7/47.6 
14.3/28.6/42.9 
11.2/22.3/33.5 



-75 
-83 
-90 
-70 
-75 
-83 
-89 



Notes. The AMBER instrument modes are Low-JHK (range 1.2— 
2.5 yum, spectral resolution R ~ 30), MR-K 2.1 (1.92-2.26 jum, R ~ 
1500), and MR-K 2.3 (2.12-2.47 /im, R ~ 1500). The baseline denotes 
the projected baseline lengths, and the position angle (PA) the projected 
position angle on sky (east of north). 



We obtained raw visibility and closure phase using the 
amdlib data reduction package dTatulli et al.ll2007t Ichelli et al.1 
2009), versi on 3.0.3. We performed a n absolute wavelength cali- 
bration as in I Wittkowski et all d201 ll) . Calibrated visibility spec- 
tra were obtained by using an average of two transfer function 
measurements taken before and after each science observation. 

Finally, we took into account that VY CMa shows an ex- 
tended dust shell that exceeds the AMBER field of view (FoV). 
In this case, the calibrated visibility values are normalized to 
the flux within the AMBER FoV and not to the larger flux of 
the more extended target, which increases the fr action of the 
correl ated flux and thus the visibility values (cf. iDriebe et al.l 
2009 ). Here, this scaling fa ctor was estimated using the results 
by Wittkowski et al. (2008) based on bispectrum speckle inter- 
ferometry with a sufficiently large FoV. They obtained a rela- 
tive flux contribution of an unresolved component of 9%/50% 
at wavelengths of 1.28^m/2.17/im. With our longer baselines, 
their extended component is fully resolved, and our visibil- 
ity data describe their unresolved component. Thus, we used a 
model of a central uniform disk (UD) and a fully resolved com- 
ponent, and scaled our visibility values such that we obtained the 
same relative flux contribution of our central component. The 
resulting scale factor for the visibility modulus was 0.45/0.90 
in the J/K bands, and we used an interpolated //-band value of 
0.65. We note that this procedure is applied in order to remove 
the instrumental signature of the limited FoV, which may be use- 
ful for a future combination of our data with data obtained by 
other instruments. However, the choice of the scaling factor does 
not affect the results presented in the following sections, except 
for the fraction of the flux that is attributed to the central stellar 
component, as the visibilities are fitted with a two-component 
model where the flux ration is a free parameter (see Sect. |4}. 



3. Atmospheric structure and morphology 

Figure [TJ shows the resulting H- and /T-band squared visibility 
amplitudes and closure phases of VY CMa obtained with the 
low resolution mode of the AMBER instrument. The /-band val- 
ues are not shown because of their lower quality . Also indicated 
are th e positions of H2O and CO bands after lLancon & Wood! 
(2000). The visibility curves show different heights, correspond- 
ing to the different projected baseline lengths. All curves ex- 
hibit a characteristic bumpy shape that resembles that of recent 
AMBER observations of Mira variables and that are interpreted 
as being indicative of the presence of mole cular layers lying 
above the continuum-forming photosphere dWittkowski et al.l 
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Fig. 1. Squared visibility amplitudes (top) and closure phases 
(bottom) of VY CMa obtained with the low resolution mode of 
the AMBER instru ment. Als o shown are the positions of CO 
and H2O bands after lLancon & Wood (2000). The different col- 
ors indicate our different ranges of PA (cf. Table [1}. 



2008, 2011). Our VY CMa visibility curves show maxima near 
2.25 yum and 1.70/mi, where the intensity of molecular bands is 
lowest. The visibility decreases toward the water bands centered 
at ~ 1.5 /mi and ~l.9fj.rn. It also exhibits a steep drop at the on- 
set of the CO absorption feature at 2.3 /mi. The decrease of the 
visibility indicates an increased contribution from extended in- 
tensity at these molecular bands. The presence of molecular lay- 
ers in the extended atmospheres of red supergiants is consistent 
with earlier observations of red supergiants such as // Cep and 
Betelgeuse dPerrin et al.l 2004, 2005; Tsuji 2006) and with the 
detection of warm water layers in the ISO and Herschel spectra 
of VY CMa dPolehampton et alfeOloHRover et alJl2010h . 

Our VY CMa closure phases exhibit complex curves that 
correlate with the positions of the molecular bands and with 
the projected position angles of the observation. Non-0/non- 
180 deg closure phase values are indicative of deviations from 
point symmetry. At the near-continuum band close to 2.25 /mi all 
curves show relatively small closure phases with absolute values 
< 20 deg. This indicates relatively small deviations from point 
symmetry of atmospheric layers close to the photosphere. These 
deviations may be caused by photospheric convection cells or 
by asymmetric molecular emission that may not be completely 
absent at a near-continuum bandpass. At the molecular bands 
of H2O and CO, the closure phase values increase, which indi- 
cates a stronger deviation from point symmetry, and also show 
variations as a function of position angle. Along the EW orien- 
tation (PA -90 deg), the closure phase values remain at abso- 
lute values below ~ 20 deg within the CO band and the water 
band. Toward the NW-SE orientation (position angle -70 deg), 
they show a ~ 90 deg signal at the center of the water band at 
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Fig. 4. Squared visibility amplitudes of VY CMa as a function 
of spatial frequency obtained at the near-continuum bandpass at 
2.20-2.25 yum. Also shown are fits of a UD (red) and a PHOENIX 
atmosphere model (blue). These contribute 50% to the total flux, 
and the remaining part is attributed to an over-resolved dust com- 
ponent. 



1.9/im, and absolute values up to ~40deg at the CO band be- 
yond 2.3 fim. The increased asymmetry of the molecular layers 
from the E-W direction toward the NW-SE direction coincides 
with the larger-scale asymmetry of the dusty circumstella r en- 
velope (PA 160deg) reported bv IWittkowski et all (Q998) at a 
wavelength of 2. 1 7 /mi. Along the PA of the stronger asymmetry, 
the sign of the closure phase differs between the H2O band CO 
band indicating that the H2O and CO regions are not co-spatial. 

Figures [2] and [3] (electronic edition) show the results ob- 
tained with the medium resolution modes at PA -89 deg and 
-83 deg. These results are consistent with the low resolution re- 
sults shown in Fig. Q] The visibility shows a maximum at the 
near-continuum bandpass at ~2. 20-2. 25 /mi and decreases to- 
ward the water band centered at 1 .9 /mi and the CO lines between 
2.3 /mi and 2.5 //m, in particular at the CO band-heads. The 
corresponding UD curve shows minimum extensions at 2.20- 
2.25 /mi and increases by about 20% at the water band near 
2.0 /mi and by up to ~50% at the CO band-heads. The closure 
phase shows relatively small absolute values of about 10 deg at 
the near-continuum bandpass, of 30^-0 deg at the water and CO 
bands, and peaks at the CO band-heads, confirming asymmetric 
structures of extended CO band-head emission. Also shown in 
Figs. [2] and [3] are visibility fits of a UD model and a PHOENIX 
model atmosphere to the 2.20-2.25 //m near-continuum band, as 
described in the following section. 

Overall, our data may q ualitatively be consistent with the 
scenario proposed by Smith et al. (2009), consisting of a cone- 
like envelope toward the SE direction that extends from a large- 
scale dusty envelope to a small-scale molecular water and CO 
envelope close to the photosphere, and additional close CO 
cloudlets that may not be co-spatial with the water envelope. 
However, for a definite and quantitative picture of the inner 
molecular envelope of VY CMa, an interferometric imaging 
campaign using a well-filled Mv-plane will be required. 



4. Fundamental parameters 

The most direct method to derive the effective temperature is 
a measurement of the radius where the Rosseland-mean optical 
depth equals 2/3 together with a measurement of the bolomet- 
ric flux. The luminosity is most directly estimated by measure- 
ments of the bolometric flux and the distance. The bolometric 
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Fig. 6. Location of VY CMa's revised fundamental parame- 
ters in the HR diagram. Also shown are evolutionary tracks 
from lEkstrom et all (120121) for masses of 15 M Q , 20 M , 25 M , 
32 M , and 40 M . The solid lines denote the models without 
rotation, and the dashed lines those with rotation. 



flux of VY CM a can be we l l estim ated using the high-precision 
photomet ry bvlSmith et al] (1200 lb complemented by the IRAS 
fluxes (Beichman et al. 1988), cf. also iJones et alj (12007) and 
IChoi et all (120081) . We de-reddened these data using an £g_y 
value of 0.6 obtained from the COBE maps (ISchlegel et al.l 
1998), and obtained an integrated bolometric flux of /t, i = 
(6.3 + / - 0.3)10~ 13 W/cm 2 . The distance estimate to VY CMa 
has recently been improved by high precision parallax mea- 
surements by IChoi et all (120081) and IZhang et all (120121) . Here 
we use the average of these two measurements, and adopt a 
distance of d - 1.17^Q^kpc. The classic distan ce used most 
often in the literature was 1.5 kpc going back to Lada & Reid 
(1978) and the suspected membership to the open cluster NGC 
2362, b ut whose distance has also been revised from 1.5 kpc to 
1.2 kpc dMel'nik & Dambisll2009h . Our adopted bolometric flux 
and distance correspond to a luminosity of L = (2.7 + QA)1Q 5 L G . 
This value is slightly smaller than recent estimates bv lChoi et al.l 
(2008) and Mauron & Josselinl (1201 ll) . mostly because of the 
slightly larger adopted distance tak i ng int o account the addi- 
tional mea surement by IZhang et a l. (2012), but otherwise well 
consistent. IChoi et a l. (2008) discussed that this value is also 
consistent with most earlier values available in the literature, if 
they are re-scaled to the same distance. We note that the error of 
L is still dominated by the uncertainty of the distance. 

The angular diamet er of VY CMa was measured 
by iMonnier et all (120041) in the K-band (2.2+0.3 /mi) to 
(18. 7+0. 3+0. 4) mas. They note that the relation of this diameter 
to a true photospheric diameter relies on addition al assumption 
on the limb-darkening, and they also mention in iMassev et al.l 
(2006) that this value may overestimate the photospheric di- 
ameter because of a possible contamination by more extended 
molecular layers and dusty features. Our spectrally resolved 
visibility data in Figs. QJ[3] confirm the presence of extended 
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layers of H2O and CO, so that a broad-band diameter estimate is 
indeed significantly contaminated by extended molecular layers 
and thus overestimates the photospheric diameter. With our high 
spectral resolution, however, we are able to estimate a near- 
continuum diameter at 2.20-2.25 yum that is largely free from 
contamination by molecular layers. With our longer baselines 
(up to 48 m compared to 38 m), our data also resolve out possi- 
ble compact dusty features to smaller scales and probe spatial 
frequencies close to the visibility null, which significantly re- 
duces calibration uncertainties. Figure 2]shows all our visibility 
data taken at spectral channels in the near-continuum bandpass 
between 2.20 //m and 2.25 yum as a function of spatial frequency, 
together with fits of a UD and a PHOENIX model atmosphere. 
The PHOENIX model was computed with v ersion 16.03.01 of the 
PHOENIX code dHauschildt & Baronll999h . We used solar metal- 
licity, mass M = 20 M Q , r eff = 3500 K, logg = -0.43. We used 
different outer pressure boundary conditions of 10~ 7 dyn/cm 2 
and 10~ 12 dyn/cm 2 , giving virtually identical results. The best-fit 
diameter using the UD model was @ud = (1 1.0 ± 0.3) mas, 
and using the PHOENIX model Ross = (11.3 + 0.3) mas (cor- 
responding to the model layer where tr oss = 2/3). Variations 
of the model mass, T e s, and log g lead to diameter variations 
of less than 0.1 mas. The flux contribution of the stellar source 
was a free parameter resulting in a value of 5 1 %. This parameter 
accounts for the fraction of an over-resolved flux component 
within our FoV and also for any residual of the uncorrected flux 
outside the FoV (Sect. [2j. The scatter of the visibility values in 
Fig. [4] may be caused by calibration uncertainties (in particu- 
lar the points at 60-65/arcsec), or by different contributions of 
an over-resolved dusty component at different position angles. 
However, the diameter remains well constrained by the mea- 
surements close to the visibility null. Fig. [5] (electronic edition) 
shows the corresponding fit using a near-continuum bandpass 
in the //-band at 1.70-1 .75 /an, resulting in a well consistent 
angular diameter of 0r oss = (11.4 + 0.3) mas, which increases 
the confidence in our diameter estimate. Figs. [2] and [3] also in- 
dicate the UD and PHOENIX models compared to the spectrally 
resolved visibility data. The PHOENIX model curve is very close 
to a UD model, and does not predict the large extensions of the 
molecular layers. PHOENIX models pr edict near-infr ared spec- 
tra of red supergiants reasonably well (Lan con et al.ll2007l and 
Figs.|2]and[3]l, which likely means that opacities from molecular 
layers are reasonably well included, but that the extension of 
the model atmosphere is currently too compact. This effect is 
most likely caused by missing pulsation and/or wind accelera- 
tion. Since our estimate of 0r oss is based on a near-continuum 
bandpass, where the contamination by molecular layers is small, 
we do not expect this to strongly affect our estimate of the pho- 
tospheric diameter. However, we can not exclude that extended 
molecular layers also contribute at our near-continuum band- 
pass to a small extent, and that our estimate of ©r oss may still 
overestimate the photospheric diameter to a small degree (of 
the order of 1 cr), so that it may represent an upper limit of the 
true photospheric diameter. The resulting effective temperature, 
based on Orqs s and fhni, is T e ^ = 3 490 ± 90 K. Using the recent 
calibration by Ivan Belle et al.l (120091) . this value corresponds to 
a spectral type of M4, w hich is consistent with the sp ectral clas- 
sification as M3/M4 bv lHouk & Smith-Moorel (Q~988) as used in 
the Simbad database. If 0r oss is considered an upper limit, the 
estimate of T e g is a lower limit. 



Table 2. Revised fundamental properties of VY CMa. 



Parameter 



Value 



Reference 



Distance (kpc) 1.17 ±0.08 1 

Bolometric flux (lO^ 13 W/cm 2 ) 6.3 ± 0.3 2 

Rosseland angular diameter (mas) 1 1.3 ± 0.3 3 

Rosseland radius (R 3 ) 1420 ± 120 4 

Effective temperature (K) 3490 ±90 5 

Spectral type M4 (M3-M4.5) 6 

Luminosity (10 5 L o ) 2.7 ± 0.4 7 

Initial mass (M ) 25 ± 10 8 

Current mass (M ) 17 ± 8 9 

Surface gravity (log g/ cgs) -0.6 ± 0.4 10 



1 We computed a large model grid for different values of M, T e g, and 
log g, which will be presented in detail in a forthcoming publication. 



Notes. 1: Averag e of IChoi et al.l d2008h and IZhang et al.l d2012h : 2 
Photometry from Smith et al. (2001) and IRAS fluxes; 3: This work 
4: 3 with 1; 5: 3 with 2; 6: 5 with Ivan Belle et all J2009T): 7: 2 w ith 1 
8-10: 5 and 7 with evolutionary tracks from Ekstrom et al. (2012). 

Figure [6] compares the resulting position of VY CMa in the 
Hertzspmng-Russel l (HR) diagr am to recent evolutionary se- 
quences by Ekstro m et al.l (|2012|) . It is located close to the red 
limit of recent evolutionary tracks of initial mass 25 M Q with 
rotation (current mass 15 M ) or 32 M Q without rotation (cur- 
rent mass 19 M ). The position of VY CMa is drawn slightly to 
the right of the Hayashi limit. This may be caused by VY CMa 
not being in full hydrostatic equilibrium. However, a slight re- 
maining underestimate of T e s, which we can not fully exclude, 
could also explain this discrepancy. Table [2] summarizes our re- 
vised fundamental properties of VY CMa. With these parame- 
ters, VY CMa would most likely not explode at its current stage 
as a red supergiant, but be at a stage before moving blueward in 
the HR diagram. 
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Fig. 2. Normalized flux, squared visibility amplitudes, UD di- 
ameters, and closure phases (from top to bottom) of VY CMa 
obtained with the MR-K 2.1 fim mode (PA -89deg). The top 
panels also show fits of a UD (red) and a PHOENIX atmosphere 
model (blue). In these model visibility fits, the stellar component 
contributes ~50% of the total flux, and the remaining part is at- 
tributed to an over-resolved dust shell. The observed and model 
flux spectra show different slopes because of the additional dust 
shell that is not present in the model, and we have subtracted 
a best-fit linear function for a better comparison of the spectral 
features. The slopes that have been compensated are indicated 
by the dashed lines. 
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Fig. 3. As Fig. E] but for the MR-K 2.3 /mi mode (PA -89 deg). 
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Fig. 5. As Fig. |U but for a near-continuum bandpass in the H- 
band at 1.70— 1.75 jum. 



